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Diffusion may obliterate fluctuation signals of the QCD phase transition in nuclear collisions at 
SPS and RHIC energies. We propose a hyperbolic diffusion equation to study the dissipation of net 
charge fluctuations. This equation is needed in a relativistic context, because the classic parabolic 
diffusion equation violates causality. We find that causality substantially limits the extent to which 
diffusion can dissipate these fluctuations. 
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I. INTRODUCTION 

Net-charge fluctuations are measured in nuclear col- 
lisions by many RHIC and SPS experiments lj. Con- 
served quantities such as net electric charge, baryon num- 
ber, and strangeness can fluctuate when measured in lim- 
ited rapidity intervals. These fluctuations occur mainly 
because the number of produced particles varies with 
each collision event due to differences in impact parame- 
ter, energy deposition, and baryon stopping. A variety of 
interesting dynamic effects can also contribute to these 
fluctuations |2|. In particular, fluctuations of mean p t , 
net charge, and baryon number may probe the hadroniza- 
tion mechanism of the quark-gluon plasma Q. 

Fluctuations of conserved quantities are perhaps the 
best probes of hadronization, because conservation laws 
limit the dissipation they suffer after hadronization has 
occurred 0, H- This dissipation occurs by diffusion. 
While the effect of diffusion on charge and baryon fluc- 
tuations has been studied in refs. 0, [g , the classic diffu- 
sion equations used are problematic in a relativistic con- 
text, because they allow signals to propagate with infinite 
speed, violating causality. 

In this paper we study the dissipation of net charge 
fluctuations in RHIC collisions using a causal diffusion 
equation. We find that causality inhibits dissipation, so 
that extraordinary fluctuations, if present, may survive 
to be detected. In sec.|n]we discuss how the classic ap- 
proach to diffusion violates causality, and present a causal 
equation that resolves this problem. The derivation we 
include facilitates our work in later sections. We general- 
ize this equation for relativistic fluids in nuclear collisions 
in sec. IIHI Related causal fluid equations for viscosit 
and heat conduction have been introduced in and 
with different heavy-ion applications in mind. 

Our causal formulation can be crucial for the de- 



scription of net-charge fluctuations, which involve rapid 
changes in the inhomogeneous collision environment [9j- 



We turn to this problem in sees. IIVIIVII In sec. IIVI wc 
discuss fluctuation measurements and their relation to 
two-particle correlation functions. In sec. we intro- 
duce techniques from [l(| to compute the effect of causal 
and classic diffusion on these correlation functions. We 
then estimate the impact of causal diffusion on fluctua- 
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FIG. 1: Density vs. position for causal (solid) and classic 
(dashed) diffusion assuming an initial delta function distribu- 
tion. The spikes represent the right and left moving diffusion 
fronts of velocity v — ±(D/Td) 



a/2 



tion measurements in sec. IV1I These three sections are 
close in spirit to work by Shuryak and Stephanov, where 
a classic diffusion model is used El. 



II. CAUSAL DIFFUSION 

To understand why the propagation speed in classic 
diffusion is essentially infinite, recall that the diffusion 
of particles through a medium is equivalent to a random 
walk in the continuum limit. The variance of the par- 
ticle's displacement increases by d 2 = (Ax 2 ) cx At in 
the time interval At between random steps. The aver- 
age propagation speed v ~ d/ At diverges in the contin- 
uum limit, where At — * with the diffusion coefficient 
D = d 2 / At held fixed. Correspondingly, a delta function 
density spike is instantaneously spread by diffusion into 
a Gaussian, with tails that extend to infinity. 

A causal alternative to the diffusion equation is the 
Telegraph equation, 



° 2 d ~ 



(1) 



where D is the diffusion coefficient and is the relax- 
ation time for diffusion 0, 0, 0] . Signals propagate at 
the finite speed v = (D/rd) 1 ^ 2 , so that a delta function 
spike spreads behind a front travelling at v, as shown 
in fig. ^ The classic diffusion equation omits the term 
cx Td- Classic diffusion describes the matter well behind 
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the front at times t ^> r d . This contrasting behavior is 
generic of hyperbolic equations such as (JIJ , which include 
a second order time derivative, compared to parabolic 
equations, like the classic first-order diffusion equation 



equ£ 
A- 



Causality concerns are not restricted to relativistic 
quarks or pions diffusing through a quark gluon plasma or 
hadron gas. They are common - and constantly debated 
- whenever time varying diffusion and heat conduction 
phenomena are discussed |Tfij . However, for most nonrel- 
ativistic systems, causality violations are minuscule, so 
that classic diffusion can be used. This need not be the 
case for relativistic fluids produced in nuclear collisions 



To motivate and understand its limitations, we first 
derive the diffusion coefficient using the Boltzmann equa- 
tion in the relaxation time approximation; see, e.g., |16| . 
Equation (JIJ can also be obtained by the moment method 
or by sum-rule arguments ^tJ- For simplicity, we fo- 
cus on electric charge transport by a single species. Gen- 
eralization to multiple species and other currents such as 
baryon number is straightforward, except for the chal- 
lenging case of color transport . We describe the evo- 
lution of the phase space distribution / using 



dj_ 

dt 



+ V p -V/=-K/-/e), 



(2) 



where v 1 is the relaxation time, v p = p/E, and E = 



{P 2 



! }V 2 . rp^g i oca j equilibrium distribution satisfies 



f e = {cxp[(E — fi) /T] ± 1} _1 where T is the temperature 
T, \x the chemical potential, and the + or — sign describes 
bosons or fermions. 

Suppose that fi differs from a local equilibrium value by 
a small sinusoidally varying perturbation 6p,, with Sn — 
(dn/dfi)Sfj,. Then / is driven from local equilibrium by 
an amount Sf(u>, k) exp{ik-x— iut}. Equation J3J) implies 



Sf(w,k) 



iv p ■ k df c 



Mw,k), (3) 



v — i(uj — k • v p ) dfj, 
plus corrections of order v . The net charge current is 



S](lj, k) = / 6f(u, k)v p dp = -ikD(w, k)«n(w, k), (4) 
where dp = d 3 p/(2ir) 3 . The diffusion coefficient is 



D(cJ,k) = - — 
3 on 



dfe 



— k • v p ) d\i 



dp, (5) 



which is a relativistic generalization of a completely stan- 
dard kinetic theory result. For k = uj = 0, we recover the 
familiar static diffusion coefficient D = ^ _1 w^ h /3, where 
the thermal velocity i? t h = 1 for massless particles. To 
obtain quantitative results from this relaxation time ap- 
proximation, we identify v~ x with the relaxation time for 
diffusion tw obtained by more sophisticated methods, see 

e.g. [H, EE El, Q 



To obtain , we omit the k dependence in , to find 

D(u,0) = D/(l-iuT d ), (6) 

where D is the static diffusion coefficient and v^ 1 = tj. 
We then write 



(f - iLUT d )j(uj, k) = -ikDniuj, k), 



to find 



d ,, . ., . 

Td^j(x,<) +j(x,t) 



-DVn(x,i), 



(7) 



(8) 



the Maxwell-Cattaneo relation Combining (JSJ with 
current conservation dn/dt + W ■ j = yields QJ. We will 
extend this argument in sec. [V] to derive eq. (|4U|) . 

Classic diffusion follows from ((HJ when the r d term 
is negligible, a result known as Fick's law. Fick's law 
violates causality because any density change instanta- 
neously causes current to flow. Including the T d term is 
the simplest way to incorporate a causal time lag for this 
current response 01- Moreover, is self consistent in 
that it saturates the /-sum rule However, while (JIJ 
is a plausible approximation, the k — > limit Q is not 
strictly justified. Possible generalizations can include the 
nonlocal equations or gradient expansions derived from 
© and (0. If leads to substantial corrections to 
classic diffusion, then such generalizations can be worth 
considering. 

In solving Q we must impose initial conditions on the 
current that respect causality. Suppose that we introduce 
a density pulse at x = at time t = 0. Microscopically, 
particles begin to stream freely away from this point. The 
current implied by JSJ is 

j(t) = 3(0)^ t/Td - [ —e-W T *DVn(a). (9) 

JO T d 

Scattering with the surrounding medium eventually es- 
tablishes a steady state in which Fick's law holds, as we 
see from (JSJ for dj/dt = 0, but this takes a time t ^> T d - 
An assumption of no initial flow j(0) = is consistent 
with our physical picture, since there is no preferred di- 
rection for the initial velocity of each particle. Note that 
an alternative choice j(0) = — DWn(0) would imply that 
the current always follows Fick's law; the corresponding 
solutions of Q would never differ appreciably from clas- 
sic diffusion. However, this choice is not causal because it 
requires that particles "know" about the medium before 
they have had any opportunity to interact with it. 



III. ION COLLISIONS 

We now extend to study the diffusion of charge 
through the relativistic fluid produced in a nuclear colli- 
sion. The fluid flows with four velocity determined by 
solving the hydrodynamic equations d^T^ u — together 
with the appropriate equation of state. We choose the 
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Landau-Lifshitz definition of in terms of momentum 
current [24(. For sufficiently high energy collisions, the 
relative concentration of net charge is small enough that 
it has no appreciable impact on u M . Correspondingly, we 
take m m to be a fixed function of x and t. 

Following [24j , we define the co- moving time derivative 
and gradient 



D T 



u^df, and 



d" - u^u v d v 



(10) 



for the metric = diag(l, — 1, — 1, — 1). In the local 
rest frame where u 11 — (1, 0, 0, 0), these quantities are the 
time derivative and —V, where V is the three-gradient. 
The total charge current in the moving fluid is j^ ot = 
nu^ 1 + j M , where the first contribution is due to flow and 
the second to diffusion. Continuity then implies 

9,K) = -a^r. (ii) 

We assume that the diffusion current satisfies 



r d D T f + f = D\7^n, 



(12) 



which reduces to JSJ in the local rest frame. 

To illustrate the effect of flow on diffusion, we consider 
longitudinal Bjorken flow, u M = (</r, 0, 0, z/t), where 
T = (^-z 2 ) 1 / 2 , 77 = (l/2)log((i + z)/(t-z)), the density 
is a function only of r and the current is — 0, 0, j z ) 
[2^|. The continuity equation i|ll|) is then 



d_ 



ld_ 

t dr 



(rn) = -d^f. 



(13) 



To evaluate the covariant Maxwell-Cattaneo relation, we 
differentiate l|12f> to find 



Tdd^d^n + d^f = -DV 2 n, 



(14) 



where V = V^V^. We then write 



d^duf) = (d ll u v ){B v j' t ) + u v B v (d tl j' t ) 



(15) 



where the second line follows from eqs. (17) and (20) of 
ref. [25|. Then O and (JTSJ) imply 



(16) 



Together, (fT3^l and (jTB|) describe causal diffusion. 

To obtain an equation analogous to Q for the ex- 
panding system, observe that the rapidity density p = 
dN/drj — A±nr, where A± is the transverse area of the 
two colliding nuclei. If one identifies spatial rapidity 77 
with the momentum-space rapidity of particles, then p is 
observable. We combine Ijl3(l and Qlfijl to find that this 
rapidity density satisfies 



Tj. 



d 2 p dp_ D d 2 p 
dr 2 dr r 2 drj 2 ' 



(17) 



where V 2 = t~ 2 d 2 / drj 2 for longitudinal expansion. Ini- 
tial conditions for p and dp/dr at the formation time t 
must be specified. In view of the causality argument sur- 
rounding @ , we assume that the initial diffusion current 
j M = 0, so that 113|) implies dp/dr = at r = r c . Note 
that the total current nu^ is initially non-zero, since 
the underlying medium is not at rest. 

In the absence of diffusion, longitudinal expansion 
leaves p fixed, as we see from <(T3|) for j = 0. Diffusion 
tends to broaden the rapidity distribution. To character- 
ize this broadening, we compute the rapidity width de- 
fined by V = ((r? — (i]} ) 2 } = N^ 1 J ri 2 pdrj, where (77) = 
and N — J pd-q. We multiply both sides of (|T7|l by rj 2 
and integrate to find 



Td 



d 2 V 
dr 2 



&V 
dr 



2D 



(18) 



Observe that classic diffusion follows from (|18|1 for Td = 
with D fixed. In that case the width increases by 



t„ V r 



classic 



(19) 



where AV = V — V(t d ). The rapidity width indeed in- 
creases, but the competition between longitudinal expan- 
sion and diffusion limits this increase to an asymptotic 
value Voo = 2D/t q . 

We now solve (|18fl for causal diffusion to obtain 



2D r ,T ° 

AV = — / F(a, X)d\, 



where a = T /rd and 



F(a, X) = a £ 



-2 e a«- 



(20) 



(21) 



We have taken dV/dr = at r = r D , as required when 
the initial dp/dr vanishes. Figure [21 compares AV/Voo 
for classic and causal diffusion. We find that the rapidity 
broadening is always slower for causal diffusion compared 
to the classic case. These solutions approach one another 
for t ^ t ; they are within 20% for r > 1.4 t for a — 
10. The classic result (O follows from (EDJ) and f ra- 
ck = T / Td — ► 00. 

We have so far assumed that Td and D are constant, 
but this may not be the case. These coefficients can vary 
with the overall density as the system expands and rar- 
efies. As an alternative extreme, suppose that diffusion 
occurs as massless charged particles elastically scatter 
with the expanding fluid of density n to t- Kinetic the- 
ory implies that t// 1 ss (eru ro i)ntot an d D ~ Td/3. If we 
take the scattering rate (<7i> re i) averaged over species and 
temperature to be constant, but assume n to t <x t x as 
determined by entropy conservation, then Td oc D oc r. 
Realistically, particle-mass effects would reduce the rate 
of growth of Td and D, as would the temperature depen- 
dence of (cr^rei). Nevertheless, the linear growth rate 
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FIG. 2: Rapidity spread vs. time for causal and classic diffu- 
sion computed using l|2()|l and I19|l respectively. Causal curves 
are for To/t^ = 0.5, 1, and 1.5. Rapidity spreads are divided 
by the asymptotic value Voo = 2D/t . 



FIG. 3: Rapidity spread vs. time for causal and classic diffu- 
sion computed using H27L 1281 and I25t respectively. Causal 
curves are for To/t^ = 1, 1.5, and 2. Rapidity spreads are 
divided by the value 2D /t = 2a/3. 



is worth considering to illustrate how rarefaction can 
change the results. 

Including this extreme effect of rarefaction implies = 
Tci(t )t/t , so that 



r d = r/a, 



(22) 



where we now fix the parameter a = T /Td(T ) at the 
initial time. Taking D — gives 



D = r/3a, 

so that the diffusion equation i fTSjl becomes 



_ d 2 V 
'~d^ 



dV 



2 

3r' 



(23) 



(24) 



As before, classic diffusion is obtained by omitting the 
second derivative term. We find 



Ay = — In—. 

3a t 



classic 



(25) 



We see that the width now increases without bound, 
which is not surprising since our time varying parame- 
ters lUllig imply that = 2D/t oc t. 

Informed by the classic case, we write the causal equa- 
tion as 



d 2 V , dV 2 



de 2 



de 



(26) 



where 9 = ]xlt/t . When a = 1, the first derivative 
vanishes, so that 



AV = (lnT/r ) 2 /3. 



For a/1, we find 
2 



AV = 



3(a- 1) 



1 



(27) 



(28) 



For a > 1, the increase of AV is always smaller 
than classic diffusion (|25|) . In this regime, we see that 



AVy(AV) c i ass i c — > a /(a — 1) as r — > oo. The difference 
of the ratio from unity in the long time limit reflects 
the fact that the rarefaction rate r _1 and diffusion rate 
t^ 1 fa v are in fixed proportion for all time, due to 12211 . 



For a < 1, i.e. for r Q < t^, eq. (|28|l implies that AV 
increases faster than in classic diffusion for r » t„. This 
behavior is an unphysical artifact of the assumption \'2'll . 
In this regime the rarefaction rate r _1 exceeds the scat- 
tering rate r~I ~ v, so that scattering cannot maintain 
local thermal equilibrium. Our hydrodynamic diffusion 
description is therefore not applicable - one must turn to 
transport theory, as discussed in j2^. One can see this 
explicitly by solving l|17|) for p(?/, r). Linear perturbation 
analysis reveals that the solutions are unstable for a < 1. 
We emphasize that for constant coefficients the solutions 
(|19fl and (|20|l are physical for all a, as generally holds for 
Td oc t z for any z < 1. 

The width computed with varying coefficients increases 
without bound, albeit slowly. This behavior distinguishes 
(|27|l and l|28l) from the constant-coefficient widths ob- 
tained from i|19f) and (|20|l . Nevertheless, the short time 
behavior shown in fig.[j|]resembles the constant-coefficient 
case in fig. [21 In both cases the classic and causal results 
converge over longer times. Note that the results in fig- EI 
are normalized to 2D/t — 2/3a using the value of D at 

To- 

To apply these results to collisions, we must specify 
the diffusion coefficient D and the relaxation time Td- 
Transport coefficients in quark gluon plasma and hadron 
matter have been been studied extensively 0, 0, |2?J 
l2lL I22I |23| . Flavor and charge diffusion estimates in a 
plasma yield D ~ 1 — 3 fm and ~ 3D ~ 3 — 9 fm [2l| . 
The hadron gas case, which is more relevant to this work, 
is complicated by the menagerie of resonances produced 
in collisions. In ref. [23|], the charge diffusion coefficient 
was computed using a hadronic transport model, yielding 
D w 2 fm and 7^6 fm. 

We argue in sec. I VII that the short time behavior in 
figs. |21 and describes diffusion following hadronization. 
Hadrons form at a rather late time, roughly r Q ~ 6 — 
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12 fm. If freeze out occurs shortly thereafter, at r/ less 
than 20 fm, then the range r/r < 3 is important. A 
hadronic relaxation time ~ 6 fm is relevant, so that 
1 < a < 2. In this case causal diffusion gives a much 
slower spread in rapidity than classic diffusion. 



IV. FLUCTUATIONS AND CORRELATIONS 

Let us now turn to net charge fluctuations and their 
dissipation. To begin, we review key features of multi- 
plicity and charge fluctuation observables. Our discus- 
sion builds on the detailed treatment in ref. [27^. which 
stresses the relation of fluctuation observables to two- 
particle correlation functions. We extend that treatment 
by identifying the net-charge correlation function 1|34|) 
that drives dynamic charge fluctuations. In the next sec- 
tion we will show how diffusion affects the evolution of 
this correlation function. 

Dynamic fluctuations are generally determined from 
the measured fluctuations by subtracting the statistical 
value expected, e.g., in equilibrium |27| . Dynamic multi- 
plicity fluctuations are characterized by 



Raa — 



(N 2 ) - (N) 2 - (N) 

W 2 : 



(29) 



where (• • •) here is the event average. This quantity is 
obtained from the multiplicity variance by subtracting its 
Poisson value (N) . Similarly, the covariance for different 
species is 



Rab 



(N a N b )-(N a )(N b ) 
(N a )(N b ) 



(30) 



which also vanishes for Poisson statistics. These quanti- 
ties depend only on the two-body correlation function 



Tab (Vl ,V2) = Pab (Vl , V2 ) - Pa ) Pb (V2 ) 



(31) 



where Pabij]!,^) = dN ab / drjidrj2 is the rapidity density 
of particle pairs for species a and b at the respective pseu- 
dorapidities r/i and 172, and p a is the single particle rapid- 
ity density. We use the same symbol for momentum and 
configuration space rapidity, because we will take these 
quantities to be equal in later sections. In [27j it is shown 
that 



-Rah = 



(Na)(N b ) 



(32) 



These quantities are robust in that they are normalized 
to minimize the effect of experimental efficiency and ac- 
ceptance. 

The STAR and CERES experiments characterize dy- 
namic net-charge fluctuations using the robust variance 



v=R A 



R — 2R^ , 



(33) 



proposed in (53; here we use v rather than the more 
conventional notation v& yn for simplicity. To establish 
the relation of v to the net charge correlation function 



q(Vi,m) = r++ + r__ - r+_ - r_+, 
we integrate over a rapidity interval to find 



(34) 



N 2 n 



qdr\\dr\2 = N{uj q - 1), 



(35) 



of 



the 

2 



net charge is Nco q 



where the variance 

((N + — A_) 2 ) — (N + — A_) 2 and the average number of 
charged particles is N = (N + + N_). Expanding (|35|> 
yields 



fi = f + R ++ + f_R— - 2f+f_R, 



(36) 



where / + = (N + )/N = 1 — / If the average numbers 

of particles and antiparticles are nearly equal, 



40. 



(37) 



HIJING simulations show that this relation is essentially 
exact for mesons at SPS and RHIC energy (we find v and 
4f2 are equal with 2% statistical error at SPS energy). We 
mention that S7 itself was suggested as an observable in 
ref. and that PHENIX measures the related quantity 
ojq = 1 + NO.. However, ft is not strictly robust. The 
next section implies that q and, consequently, fl is of 
more fundamental interest than v. However, in view of 
(j^Tjl and the qualitative aims of fluctuation studies, we 
will regard these quantities as interchangeable. 



V. CORRELATIONS AND DIFFUSION 

To understand how diffusion can dissipate dynamic 
fluctuations of the net charge and other conserved quan- 
tities, we apply the theoretical framework developed by 
Van Kampen and others ^(J- A relativistic extension of 
these techniques will enable the computation of the cor- 
relation function (|34|l as well as statistical quantities like 
v, while introducing no additional parameters. 

We start with a single charged species of conserved 
density n(x, t) that evolves by diffusion. Consider an 
ensemble of events in which n(x, t) is produced with dif- 
ferent values at each point with probability P{n(x, t)}. 
After ref. ^(ji one writes a master equation describing 
the rate of change of P in terms of transition probabilities 
to states of differing n on a discrete spatial lattice. Diffu- 
sion and flow are described as transitions in which parti- 
cles "hop" to and from neighboring points. The Fokker- 
Planck formulation in ref. |6| can be obtained from this 
master equation in the appropriate limit. Alternatively, 
one can use the master equation to obtain a partial dif- 
ferential equation for the density correlation function by 
taking moments of P. The derivation is standard and we 
do not reproduce it here |Toj . 
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If the one body density follows classic diffusion dynam- 
ics, then the density correlation function 

r(xi,x 2 ) = (nin 2 ) - (ni)(n 2 ) - 5(x x - x 2 )(ni), (38) 

satisfies the classic diffusion equation 

^-D(V; + V|))r(x 1 ,x a ) = (39) 

[To| . This density correlation function is analogous to 
rapidity-density correlation function l|31|l . If the event- 
averaged single-particle density satisfies Q , then we can 
extend this result to causal diffusion by fourier trans- 
forming H39fl and using |JBJ, as in sec. [n] The inverse 
transform yields 

(^ + ^-U(V?+V§))r(x lj x 2 ) = 0. (40) 

We remark that the delta-function term in l|38|) implies 
that the volume integral of r vanishes when particle num- 
ber fluctuations obey Poisson statistics. This term is 
derived along with i|39|l in ref. 01 j it ensures that corre- 
lations vanish in global equilibrium. 

In a nuclear collision with several charged species, only 
the net charge density n + — n~ satisfies diffusion dynam- 
ics, since chemical reactions can change one species into 
another. We define the net charge correlation function 
g(xi,x 2 ) as 

q = ((nf - n^)(r4 - n^)) - {nt - n^){nj - n^) 
-5(xi -X2)(nf +% ). (41) 

Observe that the integral 

/gdxidxa = ((N+ - N_) 2 ) - (N+ - N_) 2 

-{N++N-) (42) 

gives the net charge fluctuations minus its value for un- 
corrected Poisson statistics. This integral vanishes in 

equilibrium. We expand (|41|l to write q = r ++ + r — 

r_| — r |_, where each r ++ and r has the form 1|38|) 

and r + _ = (n+n^) - (n+)(n^). 

To obtain evolution equations for the matter produced 
in ion collisions, we incorporate longitudinal expansion 
following sec. IHII The net-charge correlation function 
q(?7i , 772, r) is given by (|34(l . where we now take 

r ab ee (p a lP b 2 ) - {pl){p\) - S ab 6( Vl - m ){pl). (43) 

This correlation function has the same form as (|41|l with 
densities n replaced by rapidity densities p = dN/dr] oc 
tit and xi : 2 replaced by spatial rapidities 771,2- We find 
that q(r]i,r]2,T) obeys 

(^ + |:-^(Vi + Vl)) g = (44) 

for V? 2 = T- 2 d 2 /drf 12 . 



To compute the effect of diffusion on net charge fluc- 
tuations in the next section, we write (|44|l in terms of 
the relative rapidity r\ T = r\\ — 772 and average rapidity 

\ d dr 2 8t t 2 dr\ 2 t 2 dvi 2 ) 

the "2" follows from the transformation to relative rapid- 
ity rj r . To compute the widths of q(rj ri r] a ,T) in relative 
or average rapidity, one multiplies l|45(l by rj 2 or r/ 2 and 
integrates over both variables. We find 

A((r? r - ( V r)) 2 ) = 2AV(r) (46) 

and 

A(( ??Q -(r /a )) 2 ) = Ay(r), (47) 

where AV(t) is calculated using l|20() (or l|28|l if a time- 
dependent Td and D is more appropriate). 

VI. SURVIVAL OF SIGNALS 

In this section we discuss how diffusion can affect our 
ability to interpret hadronization signals. As a concrete 
illustration, we suppose that collisions form a quark gluon 
plasma that hadronizes at a time r Q , producing anoma- 
lous dynamic charge fluctuations v. The hadronization 
model of Jeon and Koch Q implies that plasma produces 
a value Nu qgp « — 3. In contrast, the value for a hadronic 
resonance gas is variously estimated in the range from 
Nv « —1 j4| to —1.7 [23. Hadronic diffusion from r Q 
to a freeze out time 77 can dissipate these fluctuations, 
reducing \v\. 

We ask whether hadronic diffusion can plausibly bring 
plasma fluctuations near the hadron gas values. While 
the estimates of ref. Q can be questioned , they will 
serve here as benchmarks. In this section we emphasize 
the impact of causal diffusion on this problem, employing 
a simplified approach tailored to that aim. Phenomeno- 
logical conclusions require a more realistic model. 

Before proceeding, we emphasize that ref. Q compares 
quark gluon plasma and hadron gas fluctuations as they 
might emerge from distinct systems at equivalent chem- 
ical potential and temperature. In contrast, the fluctua- 
tions in a hadron gas that has evolved from a plasma ini- 
tial state are fixed by charge conservation. Suppose that 
hadronization in single event produces a positive charge 
in one rapidity subinterval with a compensating negative 
charge in another. The final hadronic state maintains 
that situation - freezing in the plasma fluctuations - un- 
less diffusion can redistribute the charges between those 
subintervals. 

To compute fluctuation observables following sees. IIIII 
and0 we specify the correlation function by identifying 
spatial and momentum-space rapidity at a fixed freeze 
out proper time Tf. Observe that ISR and FNAL data 
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rapidity window A freeze out x f (fm) 



FIG. 4: Rapidity dependence of dynamic charge fluctuations 
assuming diffusion has increased the width from <jq — 1 to 
Of — 2 (bottom) and ctq = 0.25 to 07 = 0.5 (top). 

|30| can be characterized as Gaussian near midrapidity. 
Moreover, these data show that charged particle correla- 
tions are functions of the relative rapidity r\ r = r\\ — 772 
with only a weak dependence on the average rapidity 
Va = + 7 ?2)/2. Near midrapidity, these data inspire 
the form 

« ^e-^ 2 -^ 2 * 2 (48) 

for S>(r. Diffusion increases the widths a and E com- 
pared to their initial values a and E G at the hadroniza- 
tion time r Q in accord with (|46[) and H47|) . For simplicity, 
we assume that E D is sufficiently large that we can ne- 
glect the time dependence of E in Ij47(l . We point out that 
(|48|l is an exact solution of the classic diffusion equation 
on an infinite rapidity interval for Gaussian initial condi- 
tions. Moreover, l|48|l is a good approximation for causal 
diffusion, provided that the rapidity region of interest 
does not appreciably exceed a or E; see the discussion in 
sec. El 

Dynamic fluctuations v are computed by integrating 
q over an interval —A/2 < 771, 772 < A/2 corresponding 
to the experimental acceptance. For E greater than A/2 
and cr, we use (|35|) and l|37|) to estimate v 

4 ,A/2 ,A/2 

iv J-A/2J-A/2 

g ,A/2 rA/2-Va 
7V JO J-A/2+r, Q 

w Nv^ erf (A/%/8ct) , (49) 

where the total number of charged particles is N w 2/?A. 
We take the rapidity density p for each charge species 
to be uniform. The quantity Nv x — 2q (27rp 2 E 2 )~ 1 / 2 is 
the value obtained for a large rapidity window. 

The dependence of the dynamic fluctuations Nv on 
the relative rapidity interval A implied by l|49ll is shown 



FIG. 5: Decrease in dynamic charge fluctuations as a function 
of freeze out time 77 for the acceptance A = 1 and <to = 0.25. 
The gray band indicates the level needed to obscure plasma 
signals. Solid and dashed curves are respectively computed 
for constant and varying coefficients using P?jl and and 
the corresponding classic eqs. 11911 and 1251 . 

in fig. 2] for several ad hoc values of a. Also shown as 
dashed curves are the ratios of these quantities, 

Nvj_ = erf(A/78a / ) 

Nv erf (A/\/8V ) ' 1 ' 

computed for a a the smaller initial width and a/ the 
larger final width. This ratio indicates how much the 
fluctuations are reduced as the width increases from a a 
to 07. fn both figures, the change in widths are chosen 
so that they nearly "hide" initial QGP fluctuations at 
the level expected in ref. 0, since Nv^/Nv^^ » 1/3 — 
1/2. Balance function measurements are consistent with 
a width a ss 0.5 |l|. ISR and FNAL experiments suggest 
cr ~ 1 in pp collisions. 

The question then becomes: are such increases in a 
plausible in a diffusion model? Equation (|46|l implies 
that cr 2 = cr 2 + 2AV(r/). The relevant "formation time" 
is the time at which hadronization occurs. This occurs 
quite late in the evolution, roughly from r D ~ 6 to f 2 fm. 
Freeze out occurs later, perhaps as late as 77 ~ 20 fm. 
The ratio t^/to is then in the range from one to two, 
where the difference between causal and classic diffusion 
is substantial, as shown in figs. [5] and [3J We take D 1=3 
2 fm and w 6 fm from ref. [23] , as discussed in sec. IIIII 

The solid curves in fig. show the decrease of dy- 
namic fluctuations computed using (|50fl . I|46|l . and 120|l . 
The dashed curves are computed using l|28|l instead of 
(|20ll . We take the value a a w 0.25 for the initial width; 
larger values imply a smaller net reduction of Nv. A 
hadronization time r Q = 9 fm is assumed, corresponding 
to a = t q I Td = 1.5 in figs. [21 and |21 The gray band indi- 
cates the level to which Nv must be reduced to hide the 
effect of plasma fluctuations. We see that classic diffu- 
sion can bring Nv to this level if the system lives as long 
as 20 fm, while causal diffusion cannot. 
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VII. SUMMARY 

In this paper we introduce a causal diffusion equation 
to describe the relativistic evolution of net charge and 
other conserved quantities in nuclear collisions. We find 
that causal limitations inhibit dissipation. To study the 
effect of this dissipation on fluctuation signals, we obtain 
a causal diffusion equation for the two-body correlation 
function. Our approach complements the treatment in 
0, but is more easily generalized to include radial flow 
and other dynamic effects for comparison to data. 

We then use these equations to provide estimates that 
show that net charge fluctuations induced by quark gluon 
plasma hadronization can plausibly survive diffusion in 
the hadronic stage. This result bolsters our optimism 
for fluctuation and correlation probes of hadronization 
and other interesting dynamics. Our aim here has been 
to emphasize the impact of causal diffusion on charge 
fluctuations. Correspondingly, we employ an idealized 
approach tailored to that aim. Phenomenological con- 
clusions require a more realistic model, to be developed 
elsewhere. 

Our results on post-hadronization diffusion owe to the 
similarity of the lifetime of the hadronic system and the 
relaxation time for diffusion r^. In fact, RHIC data sug- 
gests hadronic lifetimes that are as short or shorter than 
we assume: HBT |32j and resonance |33j measurements 
suggest a lifetime from roughly 5 to 10 fm. On the other 
hand, our estimate of rests on transport calculations in 
ref. |23j . The importance of this problem invites further 
study of charge and other transport coefficients. 

In closing, we briefly comment on the experimental sit- 
uation. Charge fluctuation measurements give roughly 



similar values from SPS to RHIC energy jlj. Further- 
more, the STAR experiment finds Nv w —1.5 at RHIC 
|28j . While this value is close to the benchmark hadron 
gas estimates 0, |28| , we emphasize that those estimates 
were constructed for an equilibrium resonance gas as- 
suming no prior plasma stage. Flow and jet quenching 
measurements provide strong indications - if not proof - 
of plasma formation 34] . As discussed earlier, Nv can be 
very different in a resonance gas formed from a hadroniz- 
ing partonic system and, indeed, should be very close to 
plasma values. 

Extraordinary charge fluctuations should be seen, but 
are not. It may be that net-charge fluctuations produced 
by plasma hadronization are closer to the benchmark 
hadronic estimates, as suggested by Bialas [221 • What 
then? Experimentally, one can turn to the centrality and 
azimuthal-angle dependence of fluctuations and, eventu- 
ally, to correlation-function measurements as more sensi- 
tive probes of net charge fluctuations. Such studies have 
already begun to yield important information P, l35| . In 
addition, one can study baryon number and isospin fluc- 
tuations, which are more closely related to the QCD order 
parameter and, therefore, to hadronization physics and 
phase transition dynamics jjj . 
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